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The topological charge is constructed for SU(3) center vortex world-surfaces composed of ele- 
mentary squares on a hypercubic lattice. In distinction to the SU(2) case investigated previously, 
it is necessary to devise a proper treatment of the color structure at vortex branchings, which arise 
in the SU(3) case, but not for SU(2). The construction is used to evaluate the topological suscepti- 
bility in the random vortex world-surface model of infrared Yang-Mills dynamics. Results for the 
topological susceptibility are reported as a function of temperature, including both the confined 
as well as the deconfined phase. 
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1. Introduction 

The random vortex world-surface model [1-8] aims to describe the gluonic dynamics of the 
strong interaction vacuum on the basis of center vortex degrees of freedom. These vortices are 
closed tubes of quantized chromomagnetic flux which spontaneously condense in the vacuum, 
giving rise to the nonperturbative phenomena which characterize the infrared sector of strong in- 
teraction physics, namely, confinement, the spontaneous breaking of chiral symmetry and the axial 
£/a(1) anomaly. 

The random vortex world-surface model, by exploring the consequences of a simplified dy- 
namics governing a reduced, infrared effective set of (vortex) degrees of freedom, complements 
studies of center vortices extracted from lattice gauge configurations [9-15]. Such studies employ 
appropriate gauge fixing and projection procedures [9, 10, 12] to determine the vortex content of 
gauge configurations; having isolated the vortices, one can investigate to what extent, e.g., con- 
finement is caused specifically by those vortices. Conversely, one can ask whether any of the 
characteristic strong interaction phenomena persist if vortices are removed from the lattice con- 
figurations. These investigations have yielded strong support for the vortex picture of the strong 
interaction vacuum; pars pro toto, results include vortex dominance of the string tension [9, 10, 12], 
the absence of topological charge and a chiral condensate when vortices are removed [11, 12], and 
a picture of the deconfining phase transition in terms of vortex percolation properties [13]. Further 
details, including an outlook on the relation of the vortex picture to other models of the strong 
interaction vacuum, can be found in the reviews [16, 17]. 

On the other hand, the simplified infrared effective description provided by the random vortex 
world-surface model permits the exploration of a wider variety of settings. For the simplest case 
of SU (2) color, studied initially, the model was shown to exhibit both a low-temperature confining 
phase as well as a high-temperature deconfined phase [1], separated by a second-order phase tran- 
sition [4]. The predictions of the model for the spatial string tension in the deconfined phase [1], 
the topological susceptibility [2] and the (quenched) chiral condensate [3] quantitatively match cor- 
responding SU (2) lattice Yang-Mills results. Building on this initial progress, the random vortex 
world-surface model was extended to the case of SU (3) color. Studies of the confinement prop- 
erties yielded a weakly first-order deconfinement phase transition [4], also seen in the vortex free 
energy [6], which represents an alternative order parameter for confinement. Furthermore, a Y law 
for the baryonic static potential was observed [5]. Spurred by recent investigations of Yang-Mills 
theories with a wider variety of gauge groups, aiming at a better understanding of possible con- 
finement mechanisms [18], the confinement properties of the random vortex world-surface model 
were subsequently also studied for SU (4) color [7] and Sp(2) color [8]. These studies showed that 
the vortex picture can accomodate such diverse color symmetries, while indicating the limitations 
of the very simple effective dynamics which had proven adequate in the SU (2) and SU (3) cases. 

After these detours into investigations of other gauge groups, the present work returns to the 
SU (3) case, widening the scope from the confinement properties studied previously, cf. above, 
to the topological properties. Specifically, the topological charge is constructed for SU (3) model 
vortex world-surfaces, where some additional complications arise compared to the SU (2) case. On 
this basis, the topological susceptibility, which determines the anomalously high mass of the r\' 
meson, is evaluated as a function of temperature. 



2 



Topological susceptibility in the SU(3) random vortex world-surface model 



Michael Engelhardt 



2. SU(3) center vortex color structure 

Center vortices are closed tubes of quantized chromomagnetic flux in three space dimensions; 
correspondingly, they are described by (thickened) closed two-dimensional world-surfaces in four 
space-time dimensions. The quantization of flux is determined by the center of the gauge group. If 
one evaluates a Wilson loop Wc corresponding to a path C encircling an SU (3) vortex, one obtains 
one of the two nontrivial center elements of SU (3), 

W c = jTr^exp (ij> A^dx^j = exp(±27n'/3) . (2.1) 



The two center elements in (|2JJ) are complex conjugates; thus, there is only one type of vortex 
flux, the two possible orientations of which generate the two signs in ( |2.1[ ). This color structure 
nevertheless admits vortex branching, in contradistinction to the SU (2) case; a vortex associated 
with a center phase exp(— 2ni/2>) = exp(4?r//3) can branch into two vortices which are each asso- 
ciated with the phase exp(2?r//3). This respects the Bianchi constraint (continuity of flux modulo 
Abelian magnetic monopoles). The presence of branchings introduces an additional complication 
for the evaluation of the topological charge of model vortex world-surfaces, an issue which deserves 
further elaboration: Evaluating the topological charge Q via 

Q = 32^2 / ^ xe iivXt Tr F mv F At (2.2) 
implies that a more specific description of the vortex gauge field must be given than just the Wilson 



loops (2.1). As a first step, it is useful to cast vortices in an Abelian gauge, i.e., the gauge field 
and the field strength F^ v are diagonal in color space. This, however, does not yet completely 
specify the color structure. One can further fix the gauge by allowing the vortex gauge field to be 
proportional only to color matrices T from a certain set. Consider the following (minimal) choice, 

re{±diag(l,l,-2)} , (2.3) 

so that there is one, and only one, generator for each of the possible vortex fluxes, 

W = (1/3) Tr exp(27nT/3) = exp(±2?n/3) . (2.4) 

Consider now what happens at a branching of a vortex associated with T = diag (— 1, — 1,2) into 
two vortices each associated with T = diag (1, 1,-2), cf. Fig. [T] (left). Flux at the branching is only 
continuous modulo 2k in each color component, since 

2k 2k 

ydiag(-l -1,2) = 2- — diag (1,1, -2) - 2rcdiag(l,l,-2) . (2.5) 

Therefore, there is necessarily a source or sink of Abelian magnetic flux, i.e., an Abelian magnetic 
monopole associated with the branching. As will become more clear further below, this intro- 
duces an ambiguity into the evaluation of the topological charge which can only be remedied by 
deforming the monopole world-line away from the branching. The minimal choice of color ma- 
trices considered above is clearly not flexible enough to allow for such a deformation. Consider, 
therefore, a non-minimal choice [19] such as 

re{±diag(l,l,-2), idiag (1,-2,1), ± diag (-2, 1, 1)} . (2.6) 
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Figure 1: Vortex color structure at branchings for a minimal choice of the set of allowed color generators T 
(left) and a non-minimal one (right); cf. main text. 



In such a gauge, one has the freedom to disassociate monopoles from branchings, cf. Fig. [l] (right); 
there, 

y diag (-1,2,-1) = ydiag(l,l,-2) + y diag (-2, 1,1) , (2.7) 

without the presence of any monopoles at the branching. This property makes the non-minimal 
description (O]) suited for an unambiguous evaluation of the vortex topological charge, as opposed 



to the minimal description (2.3). Vortex world-surfaces in the following are thus described as 



consisting of patches which are each associated with a specific color direction T from the set (2.6). 



3. Center vortex topological charge 

Vortex topological charge is generated by two geometrical features, namely, vortex world- 
surface intersection points and vortex world-surface writhe 1 [2, 19-22]. As described further in 
section 0, vortex world-surfaces will, for practical purposes, be modeled as consisting of elemen- 
tary squares on a hypercubic lattice. While for continuum surfaces, topological charge due to writhe 
is distributed continuously on the surfaces [21], on the aforementioned hypercubic model surfaces, 
topological charge due to writhe is concentrated onto lattice sites (as is, of course, topological 
charge due to vortex intersection points). Thus, at first sight, it would seem that the topological 
charge of vortex surfaces composed of elementary squares can be evaluated simply by considering 
all lattice sites in turn, at each site finding all pairs of mutually orthogonal vortex squares attached 
to the site, and weighting their contribution to the topological charge appropriately with Tr(7\r2), 
where T\ and T2 are the color directions, cf. ( |2.6| ), associated with the two squares making up the 
pair in question. However, the hypercubic description introduces two types of ambiguities which 
must be resolved before the topological charge can be evaluated in this manner. 

First, whereas continuum surfaces generically intersect at points in four dimensions, random 
hypercubic lattice surfaces share whole line segments with a finite probability, e.g., as displayed in 
Fig. |2[ Such features are artefacts due to the discrete set of directions available to the surfaces on a 



'This is easily understood by noting that, in three dimensions, a static magnetic flux in the 3-direction is associated 
with the 3-component of the magnetic field, 63, and therefore with the field strength component F12. Expanding to 
space- time, this means that a vortex world-surface running in the 3-4 directions is associated with F\2 and, conversely, 
a world-surface running in the 1-2 directions is associated with F34. If the two intersect, there is a point in space-time 
which simulta neou sly supports nonvanishing Fyi and F34, and at which topological charge density is therefore generated 
according to (2.2). Similarly, if a single vortex world-surface writhes in a way which explores all four space-time 
dimensions, topological charge density is generated. 
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Figure 2: Composing vortex world-surfaces from elementary squares on a hypercubic lattice introduces 
ambiguities, such as shared links between surfaces, which are absent in generic continuum surface configu- 
rations. 

lattice. They complicate the identification of intersection points. When whole links are shared be- 
tween surfaces, in general there is no unambiguous assignment of the attached elementary squares 
to two distinct surfaces; as a consequence, the question of whether the surfaces are actually inter- 
secting or merely touching becomes at best nonlocal, and in general ambiguous. To resolve such 
ambiguities, the hypercubic surfaces are transferred to a finer lattice, and small deformations are 
carried out until there are no lattice links left to which four or more vortex squares are attached 
(links with three attached squares of course constitute bona fide branchings). The deformations are 
such that the deformed surfaces never deviate from the original ones by more than one-half of the 
original lattice spacing; therefore, on the scale of the cutoff of the infrared effective vortex dynam- 
ics, the original and the deformed surfaces are indistinguishable. In particular, all Wilson loops de- 
fined on the original configurations remain unchanged. Nevertheless, this procedure may introduce 
additional spurious topological charge lumps on the scale of the finer lattice, and the topological 
susceptibility may therefore be overestimated as a result. The corresponding systematic downward 
uncertainty associated with the measured topological susceptibility, cf. Fig. || below, is quantified 
by scaling with the appropriate power of the vortex density, which is expected to roughly track the 
topological density, and which increases appreciably in the course of the deformation procedure. 
This issue constitutes the chief source of uncertainty in the results of this work. 

The second ambiguity present in hypercubic surface configurations was already alluded to in 
section |2[ If Abelian magnetic monopoles are rigidly tied to branchings, their world-lines neces- 
sarily run through lattice sites. On the other hand, topological charge density is concentrated onto 
the sites in the hypercubic description. This coincidence of the two also creates an ambiguity; 
at monopoles, the color direction of the vortex surface changes, and consequently, by having a 
monopole precisely intersect a site carrying topological charge density, the latter can be arbitrarily 
modified at that site. Such precise coincidence is again an artefact of the hypercubic description, 
and this ambiguity is resolved by deforming all monopole lines such that they never intersect sites 
carrying topological charge. This is possible using the vortex color structure defined by (|L6|). 

4. Random vortex world-surface model 

The assumptions underlying the random vortex world-surface model have been discussed in 
detail elsewhere [1,4]; here, a short synopsis shall suffice. Vortex world-surfaces are modeled 
by composing them of elementary squares on a hypercubic lattice. An ensemble of such surfaces 
is generated using Monte Carlo methods, respecting the (Bianchi) constraint that vortex surfaces 
be closed at each update. The spacing of the lattice is a fixed physical quantity, representing the 
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Figure 3: Fourth root of the topological susceptibility measured in the SU(3) random vortex world-surface 
model, as a function of temperature. Uncertainties due to statistical fluctuations of the random surface en- 
semble are smaller than the filled circle symbols displaying the measured values; the downward uncertainty 
shown is a systematic one, discussed in detail in section 3. 

ultraviolet cutoff of this infrared effective model. Physically, it mimics vortex thickness; e.g., two 
parallel thick vortices must remain at a minimal distance in order to still be resolved as distinct. If 
one sets the scale by using the string tension value o = (440MeV) 2 , the lattice spacing is 0.39fm. 
Finally, the vortex ensemble is governed by a world-surface curvature action, symbolically, 



i.e., each instance of two vortex squares sharing a link, but not lying in the same plane costs an 
action increment c. The lone dimensionless parameter c of the model is fixed to c = 0.21 by 
reproducing the relation between the deconfinement temperature and the zero-temperature string 
tension known from SU (3) lattice Yang-Mills theory [23], T c / -fd = 0.63. 

Generating a random vortex world-surface ensemble in this manner, and determining the topo- 
logical charge of each configuration as indicated in section || leads to the result for the topological 
susceptibility % = (Q 2 ) /V (where V denotes the four- volume) depicted in Fig. ||; the temperature 
is varied as usual by changing the extent of the Euclidean time dimension of the lattice, and the 
downward uncertainty of the measurement was estimated according to the discussion in section |3[ 

5. Conclusions 

The raw data for the topological susceptibility collected in the SU (3) random vortex world- 
surface model, cf. Fig. [3[ are somewhat higher than measurements in SU (3) lattice Yang-Mills 
theory [24] (cf. in particular Table 1 therein), which roughly indicate a value of around 190 MeV 
for the fourth root of the susceptibility at zero temperature, with values as high as 213 MeV and 
as low as 168 MeV reported; for finite temperatures, cf., e.g., [25]. However, there is a substantial 
systematic downward uncertainty inherent in the measurement of the vortex model topological sus- 
ceptibility carried out in the present work. This uncertainty is a consequence of ambiguities intro- 
duced into the world-surface configurations by the hypercubic description, cf. section ||, and of the 
procedure employed to resolve these ambiguities. Taking this uncertainty into account leaves room 



S = 



c x 



(4.1) 



6 



Topological susceptibility in the SU(3) random vortex world-surface model 



Michael Engelhardt 



for the possibility that the topological susceptibility of the SU (3) random vortex world-surface 
ensemble is consistent with SU (3) Yang-Mills theory. It should be stressed that the aforemen- 
tioned ambiguities are not intrinsic to the vortex picture, but result from the specific hypercubic 
realization of the vortex world-surfaces adopted here for practical purposes. An alternative model 
description which removes the associated uncertainty could, e.g., be achieved by constructing the 
world-surfaces as random triangulations. In this way, by letting surfaces extend into arbitrary direc- 
tions in four-dimensional space-time, the ambiguities in question are avoided. On the other hand, 
such a description is technically considerably more involved, e.g., in terms of generating a grand 
canonical ensemble of surfaces, or also in terms of the book-keeping of vortex locations. 
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